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Abstract 

The  effect  of  a  favorable  pressure  gradient  on  the  turbulent  flow  structure  in  a  Mach 
2.9  boundary  layer  {Refm  ^  1.5  x  lO*^)  is  investigated  experimentally.  Conventional  flow 
and  hot  film  measurements  of  turbulent  fluctuation  properties  have  been  made  upstream 
of  and  along  an  expansion  ramp.  Upstream  measurements  were  taken  in  a  zero  pressure 
gradient  boundary  layer  44  cm  from  the  nozzle  throat  in  a  6.35  cm  square  test  section. 
Measurements  are  obtained  in  the  boundary  layer,  above  the  expansion  ramp,  71.5  cm 
from  the  nozzle  throat.  Mean  flow  and  turbulent  flow  characteristics  are  measured  in 
all  three  dimensions.  Comparisons  are  made  between  data  obtained  using  single  and 
multiple-overheat  cross- wire  anemometry  as  well  as  conventional  mean  flow  probes.  Con¬ 
ventional  flow  measurements  were  taken  using  a  Pitot  probe  and  a  10  degree  cone  static 
probe.  Flow  visualization  was  conducted  via  imaging  techniques  (Schlieren  and  shad¬ 
owgraph  photographs).  Results  suggest  that  compressibility  effects,  as  seen  through  the 
density  fluctuations  in  the  Reynolds  shear  stress  are  roughly  10%  relative  to  the  mean 
velocity  and  are  large  relative  to  the  velocity  fluctuations.  This  is  also  observed  in  the 
total  Reynolds  shear  stress;  compressibility  accounts  for  50  —  75%  of  the  total  shear.  This 
is  particularly  true  in  the  favorable  pressure  gradient  region,  where  though  the  peak  fluc¬ 
tuation  intensities  are  diminished,  the  streamwise  component  of  the  mean  flow  is  larger, 
hence  the  contribution  of  the  compressibility  term,  is  significant  in  the  Reynolds 

shear. 
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COMPRESSIBLE  TURBULENCE  MEASUREMENTS 


IN  A  SUPERSONIC  BOUNDARY  LAYER 
INCLUDING 

FAVORABLE  PRESSURE  GRADIENT  EFFECTS 

I.  Introduction 

1.1  Motivation 

Given  the  increasing  number  of  transport  and  military  aircraft  designed  for  super¬ 
sonic  and  hypersonic  flight,  the  lack  of  available  supersonic  and  hypersonic  wind  tunnel 
facilities  has  necessitated  heavier  use  of  Computational  Fluid  Dynamics  (CFD)  as  a  means 
to  evaluate  system  performance.  Though  CFD  techniques  have  evolved  to  the  point  where 
they  provide  a  significant  amount  of  insight  into  configuration  design,  difficulties  persist  in 
accurately  modeling  the  turbulent  flow  present  in  most  high  speed  situations.  In  particu¬ 
lar,  the  accurate  resolution  of  turbulence  requires  sufficient  resolution  to  capture  motions 
on  the  order  of  the  Kolmogoroff  (8)  length  scales  for  domains  with  length  scales  on  the 
order  of  complete  aircrafts.  Simply  put,  the  solution  of  high  speed  turbulent  flows  using 
the  turbulent  form  of  the  Navier  Stokes  equations  is  beyond  the  scope  of  present  day  com¬ 
puters.  White  (18)  cites  the  example  of  a  10  cm  thick  volume  of  turbulent  flow  over  a  1.2 
m  by  2.4  m  plate.  As  the  smallest  structure  of  the  flow  is  0.04  mm,  5  trillion  mesh  points 
would  be  needed  to  resolve  this  relatively  simple  flow.  Resolution  of  the  flow  is  only  the 
first  hurdle;  there  is  also  the  number  of  operations  required  to  solve  turbulence  problems. 
For  example,  a  flow  over  a  cylinder  oi  Re o  =  lO’^,  would  require  10^^  numerical  opera¬ 
tions.  A  computer  capable  of  one  operation  per  nanosecond  would  require  320,000  years 
to  complete  this  calculation!  Credible  complete  solutions  of  the  full  set  of  Navier-Stokes 
equations  have  been  obtained  using  a  supercomputer  on  10  million  mesh  points,  however 
such  calculations  are  restricted  to  very  low  Reynolds  numbers,  i.e.  Reg  =  225,  300,  670, 
and  1410  (18).  Therefore,  these  issues  preclude  the  direct  solution  of  the  Navier-Stokes 
equations  and  dictate  that  approximate  methods  be  employed. 


1-1 


These  approximate  methods  are  typically  based  upon  the  solution  of  an  averaged 
form  of  the  governing  equations  wherein  the  dependent  variables  are  cast  in  terms  of  a 
mean  and  a  fluctuation  about  that  mean.  The  equations  are  then  solved  for  the  mean 
values  of  the  dependent  variables  (velocity,  pressure,  density,  etc.)  provided  a  “suitable” 
closure  model  representing  the  fluctuation  quantities  is  introduced.  This  then  is  the  heart 
of  the  turbulence  closure  problem.  The  selection  of  an  appropriate  model  is  based  not 
only  upon  the  form  of  the  averaging  but  also  upon  the  importance  of  accounting  for  some 
or  all  of  the  fluctuation  quantities.  The  latter  issue  has  historically  been  driven  by  the 
following  considerations:  1)  the  relative  importance  of  each  of  the  fluctuation  terms  as 
determined  by  experiments  and  2)  the  ease  of  implementation.  The  difficulty  of  describing 
and  analyzing  turbulent  flows  led  to  the  use  of  many  different  closure  methods  applied  to 
two  basic  formulations  of  the  governing  equations;  the  Reynolds  Averaged  Navier  Stokes 
(RANS)  and  the  Favre  Averaged  Navier  Stokes  (FANS).  Closure  of  these  methods  typically 
are  developed  based  upon  empirical  data  and  to  this  degree  suffer  from  a  lack  of  versatility. 

Most  current  compressible  turbulent  models  are  extended  incompressible  formula¬ 
tions  and  are  thus  limited  in  their  applicability.  Morkovin’s  hypothesis  states  that  ”the 
turbulence  structure  is  unaffected  by  compressibility  as  long  as  the  fluctuation  Mach  num¬ 
ber  is  much  less  than  unity”(ll).  A  large  quantity  of  work  has  been  based  on  this  hy¬ 
pothesis.  However,  evidence  suggests  that  these  type  of  models  perform  poorly  for  flows 
with  adverse  pressure  gradients,  shock  waves  and  favorable  pressure  gradients(9).  In  fact, 
Morkovin  (11)  showed  that  at  the  Mach  number  of  1.77,  the  density- velocity  correlations 
are  of  the  same  order  as  the  Reynolds  shear  stress.  The  many  ad  hoc  methods  employed 
have  only  helped  to  emphasize  the  difficulty  of  modeling  these  types  of  flows. 

The  current  investigation  is  aimed  at  helping  to  fill  a  void  in  the  experimental  data 
regarding  turbulent  interactions  in  high  speed  boundary  layers  having  favorable  pressure 
gradients.  There  currently  is  very  little  turbulence  data  available  in  the  literature  for  flows 
of  this  type  (favorable  pressure  gradient,  non-separated  supersonic).  All  data  that  this 
researcher  found  was  of  gross  flow  effects.  Many  studies  provided  mean  flow  data  with 
skin  friction  values,  but  failed  to  resolve  any  of  the  crucial  turbulent  components.  Lack 
of  resolution  of  this  data  limits  the  usefulness  of  the  majority  of  the  data  for  judging 
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the  relative  importance  of  the  turbulence  terms.  As  the  use  of  computational  solutions 
increases,  the  need  for  more  accurate  models  of  flow  persists.  The  particular  area  of 
interest  is  the  intermediate  Mach  region  (M  =  3  -  5)  where  the  flow  contains  large  areas 
of  expansion/shock  wave  boundary  layer  interaction.  Most  of  the  available  empirical  data 
in  this  flow  regime  is  not  of  the  ’’building  block”  variety  and  thus  is  not  appropriate  for 
turbulence  development  (6). 

The  current  study  investigates  the  turbulent  flow  structure  over  a  gradually  sloped, 
favorable  pressure  gradient  surface.  More  specifically,  several  terms  which  appear  in  the 
FANS  and  RANS  form  of  the  N-S  equations  will  be  measured  to  provide  closure  approxi¬ 
mations.  The  lack  of  quality  turbulent  data  for  supersonic  flows  has  caused  a  great  deal  of 
frustration  in  the  aerospace  community.  Some  organizations  have  gone  to  great  lengths  to 
select  the  most  valuable  data  and  find  some  criteria  to  apply  so  that  future  experiments 
would  be  of  value  in  validating  computer  models  of  turbulent  flows.  As  concerns  have  been 
raised  regarding  much  of  the  high-speed  experimental  turbulence  data  presently  available 
in  the  literature,  this  research  was  conducted  with  the  criterion  established  by  Settles  and 
Dodson  (6)  as  a  guideline.  As  such,  a  synopsis  of  the  criteria  for  experiments  suitable  for 
high-speed  turbulence  model  development  and  validation  are  given  below: 

1.  Baseline  Applicability:  All  candidate  studies  for  use  must  be  experiments  involving 
turbulent  flows  in  either  the  supersonic  or  hypersonic  Mach  number  range  (i.e.,  M  « 
3.0  or  higher). 

2.  Simplicity:  experimental  geometries  must  be  sufficiently  simple  that  they  may  be 
modeled  by  CFD  methods  “without  enormous  difficulty”. 

3.  Specific  applicability:  All  candidate  studies  passing  this  criterion  must  be  capable  of 
providing  some  useful  test  of  turbulence  modeling. 

4.  Well-defined  experimental  boundary  conditions:  All  incoming  conditions  (especially 
the  state  of  the  incoming  boundary  layer)  be  carefully  documented.  For  studies 
claiming  “two-dimensional”  flow,  data  indicating  the  extent  of  the  spanwise  flow 
variations  should  be  provided. 
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5.  Well-defined  experimental  error  hounds:  The  experimentor  must  provide  an  analysis 
of  the  accuracy  and  repeatability  of  the  data,  or  error  bars  on  the  data  themselves. 
Further,  error  bounds  on  the  data  must  be  substantiated  in  a  quantifiable  manner. 

6.  Adequate  documentation  of  data:  Data  must  be  documented  in  a  tabulated  and 
machine- readable  form. 

7.  Adequate  spatial  resolution  of  data:  Experiments  must  present  data  of  sufficiently 
high  resolution,  compared  with  the  scale  of  the  flow  in  question,  such  that  the  key 
features  of  the  flow  are  clearly  resolved. 

The  geometry  of  the  gradually  sloped  surface  is  easily  coded;  thus  it  may  be  readily 
used  for  code  validation  in  CFD.  By  experimentally  determining  turbulence  under  more 
varied  experimental  conditions,  the  understanding  of  the  nature  of  turbulent  flows  can  be 
expanded  leading  to  more  accurate  and  versatile  modeling  techniques. 

1.2  Objective 

The  objectives  of  this  thesis  are:  The  identification  of  the  important  density- velocity 
correlations  in  a  supersonic  boundary  layer  with  and  without  a  pressure  gradient.  A 
secondary  objective  is  to  provide  a  simple  geometry  test  case  for  CFD  model  development 
and  validation. 

1.3  Summary 

Chapter  two  covers  the  details  regarding  the  various  approximate  forms  of  the  govern¬ 
ing  equations,  basic  theory  pursuant  to  turbulence  modeling.  The  experimental  apparatus 
is  detailed,  including  the  wind  tunnel  geometries,  photographs  of  test  apparatus,  test  equip¬ 
ment  descriptions,  probe  designs,  data  acquisition  equipment  specifications  are  in  chapter 
three.  Chapter  four  identifies  the  flow  conditions  surveyed  and  provides  details  regarding 
flow  visualization  and  error  analysis.  Chapter  five  is  the  presentation  of  the  results  showing 
relative  importance  of  the  various  turbulent  terms.  The  effect  of  the  pressure  gradient  is 
examined.  In  chapter  six  conclusions  from  this  work  are  presented  and  recommendations 
for  further  study  outlined. 
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11.  Equation  Development 


2.1  Forms  of  the  Governing  Equations 

Computational  fluid  dynamics  (CFD)  is  playing  an  ever  increasing  role  in  modern 
supersonic  and  hypersonic  aerodynamics.  Turbulent  flows  are  characterized  by  temporal 
and  spatial  scales  that  range  over  several  orders  of  magnitude.  Direct  analytical  and  nu¬ 
merical  solutions  to  the  unsteady  N-S  equations  for  turbulent  problems  of  practical  interest 
are  highly  unlikely  in  the  foreseeable  future.  Thus,  engineers  must  rely  on  approximate 
methods  to  provide  solutions  to  practical  high  Reynolds  number  problems.  The  two  most 
prevalent  approximate  methods  are  the  Reynolds  time-averaged  and  Favre  mass-weighted- 
time  averaged  forms  of  the  N-S  equations  (RANS  and  FANS).  In  either  case,  additional 
fluctuation  cross-correlation  terms  due  to  the  nonlinear  nature  of  the  N-S  equations  ap¬ 
pear  in  the  averaged  form  of  the  equations.  Therefore,  the  averaged  systems  have  more 
unknowns  than  equations,  requiring  that  some  form  of  turbulence  modeling  be  employed 
to  close  the  system.  In  the  past,  researchers  have  often  made  use  of  ad  hoc  assumptions  to 
achieve  closure.  While  these  simplified  turbulence  models  do  serve  a  purpose  and  have  had 
a  degree  of  success,  they  are  inadequate  for  predictions  of  compressible  flow  phenomena 
such  as  shock  wave-boundary  layer  interaction  and  expansion  fans  (9).  For  example,  the 
Mach  1.77  expansion  fan/boundary  layer  interaction  data  of  Morkovin  (12)  suggests  that 
the  up'v'  term  in  the  Reynolds  shear  stress  is  of  the  same  order  as  the  typical  incompressible 
pu'v'  term. 

The  form  of  the  closure  model  depends  on  the  type  of  averaging,  FANS  or  RANS.  The 
Reynolds  averaged  equations  are  obtained  by  un-coupling  the  instantaneous  flow  properties 
into  a  time  average  mean  value  plus  a  fluctuating  turbulent  contribution,  e.g., 

u{x,  y,  z,  t)  =  u{x,  y,  z)  -f  u'{x,  y,  z,  t) 

where  u  is  the  time  averaged  velocity,  and  u'  is  the  instantaneous  turbulent  fluctuating 
component.  The  Favre-averaged  velocity  is  given  by 


P 
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and 


u  =  U{x,y,z)  +  u''(x,y,z,t) 

where  U  is  the  mass-weighted-time  averaged  mean  velocity,  and  u”  is  the  Favre-averaged 
turbulent  fluctuation.  The  form  of  the  FANS  are  very  similar  in  appearance  to  the  incom¬ 
pressible  RANS.  This  coupled  with  Morkovin’s  hypothesis  has  led  to  the  current  virtually 
universal  trend  of  adopting  the  FANS  for  high  speed  compressible  flows.  As  a  result,  prac¬ 
tically  all  compressible  turbulence  models  represent  direct  extensions  of  incompressible  for¬ 
mulations,  where  constants  are  adjusted  and  the  density  allowed  to  vary.  Such  corrections 
do  little  more  than  correlate  the  data  upon  which  the  models  are  based  (5).  The  differences 
between  the  two  averaging  methods  is  very  important.  For  example,  by  definition  the  av¬ 
erage  Reynolds  fluctuation  is  zero;  however,  the  mean  Favre  fluctuation  is  non-zero.  Care 
must  be  taken  when  comparisons  between  numerical  and  experimental  results  are  made. 
For  instance,  questions  remain  as  to  what  type  of  averaging  is  inherent  to  Pitot  pressure 
probes.  However,  it  is  clear  that  cross-wire  anemometry  responds  to  time  averaged  mass 
fluxes  and  total  temperatures.  The  present  study  uses  modern,  multiple-overheat  thermal 
anemometry  to  measure  the  compressible  Reynolds  data  in  a  Mach  2.9{Refm  =  15  X  10®) 
boundary  layer.  Detailed  three-dimensional  surveys  of  the  Reynolds  apparent  mass,  shear 
stress,  heat  flux,  and  turbulence  intensities  were  obtained.  In  addition,  conventional  mean 
flow  probes  (Pitot  and  cone-static)  and  shadowgraph/schlieren  flow  visualization  were  uti¬ 
lized.  The  data  of  the  present  study  not  only  provides  insight  into  the  nature  of  the 
high-speed  compressible  pressure  gradient  turbulence,  but  also,  based  on  the  criteria  es¬ 
tablished  by  Settles  and  Dodson  (6),  serves  as  excellent  CFD  validation  cases. 

2.2  Modeling  Issues 

The  most  prevalent  methods  of  applying  CFD  to  compressible  turbulent  flows  utilizes 
the  Favre-averaged  form  of  the  N-S  equations.  Various  forms  of  separation  of  equations 
and  closure  have  been  applied  to  simulate  turbulence  in  these  equations.  With  this  type 
of  averaging,  the  density  fluctuations  are  inherently  included  in  the  mean  Favre  velocities. 
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thus  explicit  compressibility  terms  are  absent.  For  example,  a  paper  by  Marvin  (10) 
presented  the  following  form  for  the  governing  equations: 

Continuity: 


^  d{pui) 
dt  dxi 


=  0 


Momentum:  _ 

d(pui)  d{pUiUj)  _  dp  dijij  -  pu”u'-) 

dt  dxj  dXi  dxj 

Energy: 


(2.1) 


dph  d{phUj)  _  ^  ,  // I  ~  ph"Uj)  du- 

dt  dxj  dt  ^  dxj  ^  dxj  dxj  dxj 


Turbulent  Kinetic  Energy: 


d{pk)  d(pujk) 
dt  dxj 


-pUi  Uf, 


—duk 
dxk 


d{pulk) 

dxj 


d{Uip)  du-  d^u-Tik)  d{u'l) 

oXi  oxi  oxk  dXk 


(2.2) 


Reynolds  stress  transport: 


d{ujpu':u:) 


dt 


dxi 


d{pu-  U^U-  ) 
dxA 


dxj 


dxj 
8(<P) 

dxk 


-L  -L  I  ^ 

'  ^  dxi  '  ' 


dxk 


dx  4 


(2.3) 


I  _  j.  ) 

'  dxj  '  dxj  dxj 


Eqn  2.1  has  the  same  form  as  the  typical  incompressible  Reynolds  averaged  N-S 
equations.  Numerical  analysts  using  the  Favre-averaged  form  cite  as  a  benefit  the  simpli¬ 
fication  induced  by  the  absence  of  the  fluctuating  density-generated  terms.  The  lack  of  an 
explicit  fluctuating  density  component  in  the  turbulent  shear  stresses  can  be  one  of  the 
specific  shortcomings  of  this  model  in  that  the  means  by  which  compressibility  effects  can 
be  incorporated  into  a  turbulence  model  is  less  clear. 


2,3  Turbulence  Analysis 

The  compressible,  conservative  form  of  the  RANS  equations,  written  in  Cartesian 
coordinates  are  given  by 
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Continuity: 


&p  d{puj^p'u'^)  ^ 
dt  dxj 

Momentum: 

d{pui  +  p'u'i)  d{pUiUj)  _  d{Tij  +  r^) 

dt  dXj  dxi  dXj 

Energy:  _  _  _ 

d(peo  +7^)  d{ihau^)  ^  djuifij  +  u’tIj  -  gj  -  gj) 
dt  dxj  dXj 

where  the  compressible  RANS  turbulence  terms  or  correlations  are  written  as 

mj  =  -p'u'i 


=  -pu'iU'j  -  Uip'u'j  -  Ujp'u'i  -  p'u'iu'j  (2.5) 

qj  =  ph'^u'i  +  hoP'u'i  +  Wip'h'^  +  p'Ku'i 

were  and  ho  are  the  stagnation  conditions  (ie.,  Co  =  e  +  ^UiUi)  and  ij  =  x,y  or  z.  The 
turbulent  apparent  mass  flux  (mf ),  turbulent  shear  stress  (t?),  and  turbulent  heat  flux 
(qf)  are  associated  with  the  continuity,  momentum  and  energy  equations,  respectively.  For 
incompressible  flows,  the  apparent  mass  and  the  last  three  terms  in  the  shear  stress  and 
heat  flux  are  zero. 

Note  that  the  compressible  FANS  turbulence  terms  are  expressed  as 


m'/  =  0 


t:’  =  (2.6) 

q-  =  WX 

As  previously  noted,  the  Favre  averaged  terms  are  similar  in  appearance  to  those 
of  the  incompressible  Reynolds  averaged  equations.  Bowersox  and  Schetz  (4)  have  shown 
that  cross- wire  anemometry  is  well  suited  for  measurement  of  RANS  turbulent  terms.  In 
particular,  they  have  shown  that  the  total  shear  stress  can  be  directly  measured  for  thin 
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layer  type  flows.  Further,  if  the  effects  of  the  pressure  fluctuations  on  the  hot-wire  response 
are  small,  then  the  multiple  overheat  cross-wire  results  can  be  decomposed  into  all  of  the 
terms  in  equation  (5).  This  assumption  has  been  verified  by  Bowersox  and  Schetz  (4)  for  a 
Mach  4.0  free  mixing  layer,  and  Kistler  (7)  suggested  this  is  valid  for  supersonic  boundary 
layers  up  to  Mach  4.7. 
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III.  Test  Apparatus 


8.1  Experimental  Facility 

The  present  experinaental  work  was  performed  in  the  AFIT  Supersonic  Wind  Tunnel 
Facility,  see  Figure  3.1  and  the  schematic  in  Figure  3.2  for  the  picture  of  the  AFIT 
Supersonic  Facility.  The  wind  tunnel  is  a  combination  blow-down  and  draw-down  facility. 

3.2  Air  Supply 

The  mass  flow  capacity  of  the  Atlas  Compco  brand  compressor  facility  is  approxi¬ 
mately  0.45  kg/sec.  The  system  passes  its  air  through  two  Pioneer  Air  Systems  Inc.  model 
R500A  refrigeration-type  air  dryers,  a  centrifugal  moisture  and  particle  separator  and  a 
filament-reinforced  paper  filter  before  entering  the  wind  tunnel.  The  tank  capacity  of  the 
vacuum  system  is  approximately  20  cubic  meters.  The  three  vacuum  pumps  are  manufac¬ 
tured  by  Stokes  Penwalt,  model  212-11  MicoVac  pumps.  They  are  powered  by  7.6  metric 
hp  Reliance  electrical  motors.  Figure  3.3. 

3.3  Test  Section  Design 

The  test  section  of  the  wind  tunnel  is  6.35  cm  by  6.35  cm.  The  test  section,  diffuser 
and  nozzle  are  all  constructed  from  1.905  cm  thick  aluminum  alloy.  The  test  section  side 
walls  are  constructed  of  1.905  cm  thick  optical  grade  Plexiglas,  to  allow  for  photographic 
analysis  of  the  flow.  The  overall  test  section  is  of  a  modular  design.  Each  of  the  four  sections 
are  33.02  cm  long  with  rubber  gaskets  on  one  end  joining  sections  together.  (Hex-headed 
bolts  were  used  for  the  Plexiglas  of  the  test  section,  instead  of  standard  screws,  to  prevent 
tool  slippage  and  marring  the  Plexiglas  finish).  High  vacuum  grease  was  used  on  the  gaskets 
and  sealing  surfaces.  Side  plate  bolt  holes  of  the  test  sections  were  drilled  oversized  0.8  mm. 
The  oversized  holes  allow  the  adjustment  of  the  floor  and  ceiling  to  minimize  unwanted 
upstream  shocks  due  to  misaligned  interconnection  seams.  This  also  allowed  for  some 
adjustment,  side  to  side,  of  the  wind  tunnel  walls,  which  minimized  discontinuities  in  the 
wall  seams.  Modeling  clay  was  used  as  necessary  to  block  the  remaining  seam  leakage. 
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Figure  3.1  Wind  Tunnel  and  Computers 


Figure  3.2  Wind  Tunnel  Schematic 


Figure  3.3  Wind  Tunnel  and  Computers 
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Figure  3.4  Test  Section  Profile 


Figure  3.5  Picture  of  Expansion  Ramp 

The  expansion  ramp  portion  of  the  test  section  geometry  is  a  recessed  ceiling  shown 
in  Figures  3.4  and  3.5. 

A  gradual  slope  is  used  to  provides  a  favorable  pressure  gradient  region  while  avoiding 
fiow  separation.  The  approximate  slope  of  the  design  is  a  one  degree  change  in  slope  every 
0.8467  cm.  The  test  section  contour  matches  the  curve  fit: 

Y  =  2.54  *  («„  +  ai(a;/2.54)  +  a^{xl2.bYf  +  az{xl2.bAf)  (3.1) 

where  Oq  =  -0.0818011,  ai  =  0.0897045,02  =  -0.0240686,  and  =  -0.000232096.  This 
curve  fit  has  a  standard  deviation  of  0.00081  cm  and  x  and  Y  are  in  cm.  The  x  coordinate 
is  measured  from  the  beginning  of  the  contour. 
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3.4  Data  Acquisition 

Instrumentation  for  the  wind  tunnel  is  attached  to  a  Multipro  data  acquisition  sys¬ 
tem  manufactured  by  Nicolet  Instrument  Corp.  The  Multipro  has  four  model  120  data 
acquisition  boards  installed  with  1  megabyte  of  memory  per  board.  Bus  triggering  was 
used  on  the  Multipro  to  assure  data  acquisition  synchronization. 

3.5  Traversing  Equipment 

The  stepper  motor  used  is  a  MD-2  with  400  steps  per  revolution.  The  probe  traverse 
slide  is  made  by  UniSlide  and  has  a  total  travel  distance  of  16.51  cm.  The  slide  moves 
on  a  shaft  of  40  threads  per  2.54  cm.  The  Figure  3.6  shows  the  schematic  of  the  unit, 
S  is  6.35  cm  and  L  is  22.86  cm.  A  picture  of  the  traversing  assembly  is  provided  in 
Figures  3.7  and  3.8.  The  hand  cranks  in  the  Figures  provide  manual  control  of  horizontal 
and  vertical  alignment.  The  maximum  speed  of  the  MD-2  is  approximately  0.25cm/sec. 
The  circular  plate  is  used  to  vary  the  angle  of  the  probe  traversing  the  flow  so  that  it 
remains  perpendicular  to  the  surface  of  interest.  The  LVDT  (Linear  Voltage  Displacement 
Transducer)  used  is  TRANS  TEK  model  0217-0000  F-90  (Not  Shown).  The  power  supply 
for  the  LVDT  is  a  Hewlett  Packard  6205C.  Traverse  control  was  accomplished  using  a 


3-4 


Figure  3.7  Traverse  Assembly 


Figure  3.8  Traverse  Assembly 
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Figure  3.9  Schlieren  Camera  and  Mirror  Setups 


Zenith  386  computer.  Linking  of  the  traverse  equipment  to  trigger  the  data  acquisition 
allows  accurate  location  versus  data  tracking. 

3,6  Shadowgraph  and  Schlieren  Optics 

The  shadowgraphs  were  acquired  using  an  arc  light  source  made  by  Cordin  model 
5401  with  a  spark  duration  of  600  nanoseconds.  The  light  is  reflected  from  a  101.6  cm 
focal  length  mirror  and  photographed  with  series  52  Polaroid  film.  The  schlierens  used  the 
same  light  source  but  used  an  additional  mirror  and  knife  edge  to  split  the  incoming  light 
to  the  film.  The  diagrams  for  the  schlieren  layout  are  provided  in  Figure  3.9c. 
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5.7  Pressure  Probes 


Measurement  of  the  P^oo  Pjco  S'l’e  accomplished  using  Pitot  tube  and  thermo¬ 
couple  located  upstream  of  the  nozzle  in  the  settling  chamber.  Mean  flow  in  the  test 
section  is  acquired  by  the  use  of  a  Pitot  tube  to  obtain  the  total  pressure  (P12),  and  conical 
probe  for  the  static  pressure  (Pc)-  For  most  of  the  test  runs  signal  Altering  for  the  pres¬ 
sure  transducers  were  passed  through  an  ENDEVCO  model  4225  signal  conditioner  before 
passing  into  the  Multipro.  During  the  multiple  overheat  runs  the  pressure  transducers  are 
connected  to  an  ENDEVCO  Model  4428A  signal  processor  and  calibrator.  This  unit  uses 
the  maximum  range  of  the  pressure  transducer  and  spreads  the  range  of  response  linearly 
over  5  volts.  The  unit  is  self-zeroing  and  connects  directly  to  the  Multipro.  The  pressure 
transducers  documentation  assures  that  accuracy  is  below  a  maximum  error  of  0.5%  of  full 
scale  and  guarantees  linear  behavior.  A  1.0206  atm  transducer  is  used  for  the  Pitot  and 
cone  static  probes  with  an  ±  0.0051  atm  error.  The  upstream  transducer,  6.804  atm,  has  a 
error  ±  0.017  atm.  Atmospheric  pressure  is  recorded  using  a  Tranamerica  Delaval  digital 
barometer  model  number  2500-0101. 

The  conventional  probes,  Pitot  and  cone-static,  are  both  machined  from  stainless 
steel  and  are  curved  to  the  same  approximate  length  as  the  hot  wire  probes.  The  length 
is  matched  so  that  the  data  taken  would  be  as  close  as  possible  to  the  same  streamwise 
location.  The  original  0.64  cm  diameter  probe  holder  created  excessive  blockage  causing 
the  tunnel  to  unstart.  Figure  3.10.  This  probe  effectively  blocked  off  more  than  8.0%  of 
the  flow.  This  lead  to  the  development  of  a  variable  diffuser  and  a  smaller  probe  holder 
(0.3175  cm  diameter).  The  smaller  size  was  still  rigid  enough  to  minimize  deflections  of 
the  probe  yet  not  choke  the  flow. 

3.8  Temperature  Measurement 

The  temperature  transducer  used  is  made  by  Omega  Engineering,  Incorporated  and 
is  part  number  SEFE-K-5  and  has  an  accuracy  of  ±1°  K  (15).  The  transducer  is  connected 
to  a  temperature  display,  also  made  by  Omega,  the  Digicator  model. 
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Figure  3.10  Unstart  Condition 


Figure  3.11  UV  Hot  Film  Probe,  Model  1243-20 
3.9  Hot  Film  Probes 

The  cross  wire  probes  used  are  TSI  brand,  model  1243-20  and  1243AN-20  with  a 
platinum  hot  film  (Figures  3.11  and  3.12).  The  coefficient  of  resistance  is  0.24%/'’C’. 
They  are  connected  to  a  TSI  brand  IFA  Intelligent  Flow  Analyzer  100  with  internal  bus 
resistance  of  50  ohms.  A  Tektronix  2454B  oscilloscope  is  used  to  view  and  tune  the 
frequency  response  shape  of  the  probe.  A  list  of  resistances  of  the  probes  is  given  in 
Table  3.1,  where  (R)  stands  for  repaired  probe. 

For  multiple  overheat  work  an  external  resistor  setup  is  used  which  is  similar  to 
the  one  outlined  by  Bowersox  (2).  The  LVDT  power  light  is  wired  to  trigger  a  Wavetek 
model  278  function  signal  generator.  This  in  turn  set  the  cycling  frequency  of  the  resistors 
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Figure  3.12  UV  Hot  Film  Probe,  Model  1243AN-20 


Table  3.1  Hot  Film  Probe  Operating  Resistance 


Probe  type 

S/N 

Channel  1 

Channel  2 

1243AN 

944012 

4.80 

4.78 

1243AN 

944010 

4.80 

4.73 

1243AN(R) 

944012 

5.608 

5.031 

1243 

943011 

5.19 

5.21 

1243 

943009 

6.48 

6.25 

1243 

943008 

6.44 

6.34 

at  4  Hz.  Single  overheat  (SOH)  measurements  were  obtained  from  a  1:1  bridge  using 
an  external  resistor  to  set  the  overheat  level.  A  resistor  of  10  ohms  was  used.  When 
tested  for  accuracy,  the  two  resistors  used  measured  as  10.000  ohms.  Cable  resistance  was 
measured  at  0.4  ohms.  Thus,  the  total  overheat  was  based  on  10.4  ohm  total  resistance. 
This  level  of  resistance  was  above  that  listed  in  any  of  the  probe  literature.  The  higher 
overheat  resistance  gives  a  more  accurate  response  to  mass  flux  and  thus  less  susceptible 
to  temperature  fluctuations.  TSI  recommended  overheat  is  estimated  at  approximately 
1.6.  The  SOH  for  this  investigation,  range  between  1.6  and  2.2.  The  external  resistance 
used  for  the  multiple  overheat  (MOH)  is  listed  in  Table  3.2,  all  values  in  Tables  3.2  and 
3.1  must  have  0.4  Ohms  added  to  them  to  get  accurate  overheat  ratios. 

The  tuning  of  the  frequency  response  is  the  most  difficult  part  of  hot  wire  usage.  The 
probes  have  to  be  tuned  at  or  near  test  conditions.  The  only  available  test  source  is  the 
tunnel  itself,  which  only  has  a  usable  run  time  of  20  to  25  seconds.  Tuning  each  channel 
took  about  one  day.  The  consistency  of  the  channel  response  was  verified  multiple  times. 
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Table  3.2  MOH  Operating  Resistance 


Overheat  No. 

Resistance 

1 

10.000 

2 

9.524 

3 

9.091 

4 

8.750 

5 

8.333 

6 

7.778 

7 

7.500 

8 

7.143 

especially  when  switching  from  the  UV  to  the  UW  probes.  Channel  1  bridge  compensation 
was  set  to  1340  and  channel  2  was  set  to  740. 

During  single  overheat  runs  the  data  acquisition  is  triggered  by  the  sharp  increase  in 
tunnel  pressure  at  start.  The  pressurization  of  the  wind  tunnel  is  the  last  event  before  the 
IFA  is  started  and  the  traverse  is  started  across  the  flow.  Using  this  single  event  to  trigger 
the  other  channels  for  data  collection,  the  LVDT,  upstream  pressure,  and  two  channels  for 
the  hot  wire,  all  are  simultaneous.  Data  was  manipulated  in  groups  of  1024  points  for  the 
SOH.  Details  regarding  the  operating  procedures  for  the  facility  are  contained  in  Appendix 
J.  MOH  anemometry  is  more  fatiguing  on  the  hot  wire  probes.  The  higher  overheat  ratios 
and  the  cycling  of  overheat  ratios  are  more  detrimental  to  the  survival  of  the  hot  wire.  The 
short  run  time  of  the  tunnel  necessitates  a  very  quick  cycling  of  the  overheat  ratios.  The 
circuit  outlined  in  Bowersox  (2)  is  used  to  cycle  the  overheat  ratios.  Samples  are  taken 
while  the  traverse  pauses  for  three  seconds  in  even  increments  as  it  traveled  through  the 
boundary  layer.  The  trigger  for  data  acquisition  in  this  case  is  the  traverse  stopping.  This 
event  triggers  the  the  Wavetek  signal  generator  which  cycled  the  resistors  at  a  4  Hz  cycle 
speed  through  the  multiple  overheat  circuit.  The  sampling  rate  of  the  Multipro  is  set  at 
10  kHz.  As  a  result,  it  captures  2500  data  points  for  each  overheat.  Of  these  data  points, 
approximately  2056  centered  over  the  sample  time  were  used.  The  remaining  points  were 
neglected  to  eliminate  any  transients  due  to  scanning. 
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3.10  Computers  Used 

Data  is  stored  on  a  Compaq  386/25m,  and  a  Zenith  486.  The  traverse  equipment  is 
run  from  a  Zenith  386.  The  data  reduction  is  done  on  an  Amiga  2000  and  on  Sun  Sparc 
20 ’s.  The  Amiga  is  used  for  the  conventional  probe  mean  flow  pressure  data  breakdown. 
The  Sparc  20’s  are  used  for  all  other  data  analysis. 

3.11  Settings 

A  pressure  of  approximately  1.5  atm  (to  achieve  two  atmospheres  in  the  wind  tunnel 
due  to  back  pressure  losses)  is  used  as  an  air  source  and  is  exhausted  into  an  evacuated 
tank,  initially  drawn  down  to  approximately  18  mm  of  mercury  absolute.  The  vacuum 
system  initially  takes  15  minutes  to  draw  down,  and  then  subsequent  runs  take  from  8  to 
10  minutes  to  draw  down.  The  compressor  facility  had  no  problems  keeping  up  with  the 
demand. 

The  traverse  equipment  speed  has  a  range  of  speeds  from  0  to  32000,  0  being  the 
fastest.  Speeds  faster  than  370  are  not  used  due  to  system  response  losses.  For  the 
conventional  probe  usage  the  speed  is  set  at  1200,  which  corresponds  to  .127  cm/sec.  As 
the  hot  wire  data  is  sampled  at  a  higher  rate,  the  traversed  speed  is  increased  to  370, 
leading  to  a  speed  of  .254  cm/sec.  The  multiple  overheat  runs  are  restricted  to  the  first 
1.27  cm  of  the  boundary  layer  due  to  the  need  to  stop  and  sample  at  each  location  for 
three  seconds.  Samples  are  taken  every  0.254  cm  for  the  multiple  overheat  data. 

The  sampling  rate  of  the  Nicolet  Multipro  is  set  to  200  Hz  for  the  sampling  of  the 
conventional  probes.  The  sampling  rate  for  the  hot-wire  probes  is  set  to  10  kHz. 

The  frequency  of  the  sampling  rate  was  also  a  major  concern.  The  frequency  of  10 
kHz  was  chosen  to  capture  the  overall  energy  of  the  flow.  The  signal  was  not  filtered  and  a 
spectral  analysis  was  not  necessary  for  the  current  study  and  the  energy  of  the  flow  could 
be  captured  without  resolving  where  the  frequency  of  the  energy  originated.  Based  on  the 
work  of  Smits  and  Muck  (14),  the  energy  spectra  range  of  our  sampling  should  be  able  to 
capture  the  flow  character  very  well. 
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3A2  Error  Analysis 


Based  on  the  analysis  presented  in  (2)  a  Euclidean  (or  L2)  norm  are  presented  as 
measures  of  the  error.  The  L2  norm  is  given  by 


(3.2) 


where  i  indexes  the  various  errors  associated  with  measurement  x. 


3.12.1  Conventional  Probes.  Volluz  (17)  reports  that  turbulence  induces  about 
±0.0068  atm  error  for  both  probe  types,  Pitot  and  cone  static.  The  calibration  error 
gp,,  =0.005  atm  is  for  both  probes.  Upstream  pressure  has  an  error  of  gp,,  =0.017  atm.  The 
temperature  error  is  gT,,  =  1°  K 

The  probe  location  errors  are  estimated  as  ±0.5%.  Therefore,  the  error  in  the  Mach 
number  at  M  =  2.9  is  approximately  ±0.1  or  ±3%,  respectively.  The  error  in  the  local 
mean  mass  flux  is  ±2.0%. 


3.12.2  Hot  Wire  Probes.  The  linearizations  errors  have  been  shown  by  Bowersox 
(2)  to  be  negligible  compared  to  the  above  mentioned  calibration  errors.  Perturbation 
analysis  also  shows  that  probe  misalignment  errors  are  very  small,  primarily  due  to  the 
procedure  of  re-calibrating  the  probes  at  each  new  station  and  thereby  accounting  for  any 
misalignment.  The  following  errors  were  documented  for  this  equation  model  in  Table  3.3. 
Where  gp^  is  the  mass  flux  error,  gT<  is  the  temperature  error,  g/  is  the  f  sensitivity  error. 


Table  3.3  Approximate  System  Error 


Identifier 

Percentage 

^  pu 

2.0 

^Tt 

0.3 

e/ 

1.0 

0.3 

«()'()' 

8.33 

Cti 

4.17 

Cg  is  the  g  sensitivity  error,  cqiqi  errors  for  the  cross-wire  shear  and  total  temperature 
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flux  correlation  terms,  and  Cti  is  the  error  of  the  turbulence  intensities.  Based  upon  the 
analysis  presented  in  reference  (2).  This  results  in  a  total  hot  wire  error  of  8.5%. 
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IV,  Data  Reduction 


4  A  Use  of  Shadowgraphs  and  Schlierens 

Shadowgraphs  and  schlieren  photographs  were  taken  to  afford  examination  of  some 
general  flow  characteristics.  From  these  photographs  the  boundary  layer  thickness  was 
estimated  to  be  1.25  cm  thick  at  x  =  59  cm  (before  the  expansion).  The  thickness  of 
the  boundary  layer  over  the  expansion  ramp  was  estimated  as  1.9  cm  at  x  =  75  cm  (at 
the  end  of  expansion  ramp).  No  separation  was  observed,  thus  all  probes  use  were  angled 
(i.e.  perpendicular  to  the  surface)  assuming  that  the  flow  followed  the  contour  of  the  wall 
within  the  boundary  layer. 

4^2  Mean  Flow  Data  Reduction 

A  least  squares  curve  fit  for  the  solution  of  the  Mach  number  from  the  ratio  of  the 
pressure  obtained  from  a  10°  axisymmetric  cone  and  a  Pitot  probe  (2)  is  given  by: 

—  =  -0.052976  +  4.6840a;  -  18.678a;^  +  50.7006a;^  -  54.1577a:^  (4.1) 

M 

where  x  represents  PclPt2  to  calculate  the  Mach  number  and  is  good  in  the  Mach  number 
range  from  1.5  to  4.4.  The  curve  fit  has  a  standard  deviation  of  0.0006.  The  remaining  flow 
properties  were  computed  from  the  typical  compressible  gas  dynamic  relations  (1).  The 
data  from  the  conventional  probes  is  processed  primarily  on  an  Amiga  2000  running  at  28 
MHz.  Data  points  were  averaged  together  in  groups  of  100.  Data  obtained  using  the  Pitot 
probe  is  obtained  closer  to  the  surface  than  that  garnered  using  the  cone-static  probe.  The 
cone  static  probe  data  is  extrapolated  using  a  third  degree  polynomial  curve  fit.  To  match 
the  data,  procedures  for  interpolation  and  extrapolation  are  used  from  ’’Numerical  Recipes 
in  C”(13).  This  process  adds  two  points  to  the  Mach  number  profile  near  the  wall.  The 
lowest  data  point  for  the  hot  wire  data  is  lower  than  the  lowest  cone  static  probe  data. 
The  extrapolation  of  the  cone  static  data  below  the  cross  wire  data  proved  invaluable. 
Approximately  10-12  data  points  in  the  hot  film  data  were  added  by  this  small  addition 
of  the  first  of  these  two  points  in  the  conventional  data  (the  lowest  point  is  not  used  for 
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the  hot  wire  data  analysis).  The  profiles  for  the  upstream  test  location  were  done  with 
multiple  passes  to  insure  repeatability  from  one  run  to  the  next. 

4.3  Hot-Wire  Data  Reduction 

Hot  film  measurements  are  obtained  at  each  station.  Because  of  the  fragility  of 
the  probes,  the  second  run  at  a  given  station  is  usually  the  calibration  run.  All  hot  wire 
measurements,  UV  and  UW,  are  taken  for  a  set  configuration  assembly  to  remove  anomalies 
induced  by  configuration  changes.  Single  overheat  (SOH)  and  multiple  overheat  (MOH) 
measurements  were  obtained.  Eight  overheat  ratios  were  employed  at  each  MOH  point,  but 
in  most  cases  the  six  higher  ratios  are  used  to  reduce  the  data.  Lower  overheat  ratios  are 
less  sensitive  to  mass  flux  fluctuations  and  will  sometimes  become  oscillatory,  rendering 
their  response  meaningless.  The  oscillatory  behavior  increases  as  the  wire  temperature 
approaches  the  flow  total  temperature.  The  multiple  overheat  cross-wire  data  reduction 
techniques  developed  in  Bowersox  (2)  and  Bowersox  and  Schetz  (4)  were  used  in  the  present 
study.  A  summary  of  their  basic  formulation  is  presented  in  this  Section.  Multiple  overheat 
cross-wire  anemometry  in  supersonic  flow  provides  the  following  “conservative”  variable 
turbulence  data 

{pUiYipUj)',  T/^  {pUiyTl 

The  basic  form  of  a  hot  wire  is  a  cylinder  and  the  form  of  the  Nusselt  number  for  com¬ 
pressible  flow  in  dimensionless  heat  transfer  is: 

Nu  =  F{Lld,M,Pr,Ree,T)  (4.2) 

L/d  is  the  wire  aspect  ratio;  M  is  the  Mach  number;  Pr  is  the  Prandtl  Number;  Rce  is 
the  effective  cooling  Reynolds  number  (base  on  wire  diameter);  and  r  is  the  temperature 
loading  factor.  The  temperature  loading  factor  can  be  expressed  as  r  =  (T^,  -Tg)/?*  where 
Ty,  is  the  wire  temperature  are  the  temperature  the  unheated  wire  would  attain  if  placed 
in  the  flow,  called  the  equilibrium  temperature.  For  Reynolds  numbers  greater  than  about 
20,  Te  is  about  97%  of  Tf  When  the  Mach  number  normal  to  the  wire  is  greater  than 
1.2,  or  Msin<j)  >1;  Pr  is  constant;  and  the  aspect  ratio  >  1,  the  function  for  the  Nusselt 
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number  becomes: 

Nu  =  f{Re,,T)  (4.3) 

Following  the  assumptions  and  development  of  Bowersox  and  Schetz  (4) 

Nu  =  a\/  Rcg  +  b  (4.4) 


and  equal  to 


Nu  = 


TrktL{T^  -  Tt) 


(4.5) 


where  q, 


w  — 


'■R^  and  =  V^I{R^  +  Rs  +  Rl)  from  anemometer  circuit  analysis.  As¬ 


suming  Te  =  Tt  this  equation  reduces  to 


Nu 


V  ^R 


{Rw  +  jRs  d"  RlY  '^^tLiTyj  —  Tt) 


(4.6) 


The  power  laws  for  viscosity  and  thermal  conductivity  were  used 


kt  =  ko[Yj  =  k-o 


(4.7) 


with  Up  =  0.77  and  rik  =  0.89  (2).  Combining  the  previous  equations,  the  hot  wire  response 
equation  is 

~  =  (^y  [a/R^(rjr„)-”'‘/"  +  b]  (T^  -  Tt)  (4.8) 

Reo^  is  the  effective  Reynolds  number  with  }i  ^  and 


Co  =  +  + 

Rw 


(4.9) 


Replacing  T^,,  ReOe,  and  T  by  their  mean  plus  fluctuating  components,  using  the  binomial 
theorem,  retaining  only  the  first  order  terms,  and  noting  that 

(T„-I))  (4.10) 
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then  solving  for  V^/V^,  the  hot  wire  fluctuation  equation  is  derived  in  the  form 


K  V  ) 


(4.11) 


and  the  hot  wire  sensitivities  are 


a\  Re f 


and  g  = 


^Tt  U); 

2(T„-7^)^T 


fn^, 


(4.12) 


Rewriting  equation  4.8 


^ Re^  “I”  ^iRt\J R^e  ~i~  Vi'd^t 


(4.13) 


where  x,  =  yi  =  —biKaiT^i),  and  Zi  =  l{CiaiT^i)  —  bi/ai.  C  is  and  ko  is 

kt-  The  subscript  i  indexes  over  the  overheat  ratios.  At  a  minimum  two  wire  overheats 


are  needed  to  solve  for  the  unknowns  \J Re,,  and  Tj.  If  more  overheats  are  available  a  least 
squares  analysis  is  used  (2). 


NjRe,  +  Tt{Yl  Vi  -  E  +  ‘dTt'J Re^  E  Xi  + 


T/Z^iVi  +  Tt  = 


(4.14) 


-J:xiZi)  +  T,Zy!  +  2T,^ 

+TtJ Re e  E  =  E  Vi^i 


XiVi  + 


J2xi 


(4.15) 


N  represents  the  number  of  overheat  ratios.  Summations  cover  the  i  overheat  ratios.  The 
solution  was  reached  iteratively  using  the  secant  method.  Squaring  and  averaging  equation 
(4.11)  yields 


2/  Reoj' 


(4.16) 


Three  overheat  ratios  are  required  to  solve  this  equation.  If  more  overheat  ratios  are 
available  a  least  squares  approach  is  employed.  The  method  used  for  solution  of  this 
problem  is  Quadratic  Least  Squares  (QLS).  It  assumes  that  the  errors  in  the  hot-wires  are 
purely  random  and  all  bias  resides  in  f  and  g.  Therefore,  the  new  wire  temperature  loading 
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factor  is 


(4.17) 


(r.  -  T,) 


and  the  new  sensitivities  can  be  expressed  as 


f  =  fo  +  mj- 


g  =  go  +  rUg  — 


(4.18) 


where  fo  and  w,j  are  determined  from  wire  calibration  and  go  —  Tiki‘s  '^g 

-(.5  +  Tig^mf). 

Through  extensive  manipulation  of  variable  and  constants,  a  three  by  three  system 


of  equations  results 


(1. 

k)  =  “•  +  ”■  (;^  j  +  “Tr- 


(4.19) 


2fomj  2{fomg+gomf)  2goi 


2mfma 


and  the  temperature  is  calculated  using 


(  Reo^Y 
\  / 
(Reo^rTZ] 

V  Ti  J 


(4.20) 


+  - ^{Rn,  -  Rref)  +  273.15 

Kref  ^ref 


(4.21) 


The  response  of  the  cross  wires  are  independent  of  Mach  number  as  long  as  M  sm{(f))  >  1, 
but  they  are  still  dependent  on  the  normal  component  of  the  Reynolds  number.  To  derive 
the  formulas  for  analysis,  the  effective  Reynolds  number  must  be  used  in  the  x  and  y 
coordinate  system  and  related  to  oblique  shock  theory. 


cos(</>)  sin(<^)  /  Rex 
—  sin((?!))  cos(</))  I  Rey 


(4.22) 
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(4.23) 


Re,^  =  Ren  +  k.^Ret"^  =  A^Re^'^  +  2A‘2Re^Rey  +  Agiie/ 


where  A,  are  given  by 

Ai  =  cos'^{(j))  +  sin^((?!>) 

A2  =  (1  -  A:/)  cos{4>)  sin(^)  (4.24) 

As  =  kc'^  cos^(</>)  +  sin^(^) 


and  (j)  is  the  incidence  angle  of  the  hot  wire  to  the  flow.  Using  first  order  assumptions  k^ 
can  be  shown  to  be  (2) 

kc  =  '^k  (4.25) 

P2 


where  k  =  f(L/d) 

Ren  =  -Ren  (4-26) 

Pi 

the  subscript  “1”  refers  to  the  condition  just  prior  to  the  oblique  shock,  “2”  refers  to  the 
condition  after  the  shock  and  k  is  to  be  determined  experimentally  from  the  calibration 
data. 


Replacing  Reo^,  Reo^  and  ReOy  by  their  mean  fluctuating  components  and  applying 
the  binomial  theorem  and  using 


Ro  = 


B2  = 


Ai 

B3 

Ai 

B3 


B3  =  Ai  +  2A2R0 


(4.27) 


ReOej  =  Reo^y/B^ 


(4.28) 


<  1  is  assumed.  J  indexes  as  either  the  one  or  the  two  wire  on  the  cross-wire  probe. 
Solving  this  set  of  equations  and  decomposing  into  x  and  y  components 


_  ReOel^  fA2l—ReOc2  M22 

AII/A21—A12/A22 
1  Reoci^ fAii  —  ReOa2  IA12 
2ReOx  A2\I  A\,\—  A22I  A\2 


(4.29) 
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'  ReoJ 


ReOx 


D. 


1  /  ReoJ 


-621  \  R^Of. 


2  ^  ReoJ^  ^  ReoJ 

B21B22  \  Reoe  j  ^  \  Reoe 


+ 


B22 


ReoJ 

ReOf 


/2j 
(4.30) 


ReOy 


Reo^ 


DJ 


1  (Reoj\ _ 

BnB 


ReoJ  ReOy\  _  1 

ReOx  ReO;^  j  2jBii52i 

where  Di  {B21IB11  -  B22IB12)  and  D2  =  {BnfB2i  -  5i2/522)-  The  covariance  of  the 
two  wires  can  be  expressed  as 


ReoJ 

Reo. 


ReoJ 

ReOr 


+  Bl, 


ReoJ 

Reo^ 


(4.32) 


ReoJ 

ReOg 


+6^152 


TJ 

Jt, 


+  /lff2 


ReoJ  Tj 


ReOe  Tt 


+/251 


ReoJ  Tj 

ReOe  Tt  ^ 

(4.33) 


The  fluctuation  of  the  Reynolds  number  can  now  be  computed.  The  Reynolds  number 


total  temperature  correlations  result  in 


ReoJ  Tj  _  1  1  lReojTj\  1  f  ReoJ  Tj\ 

ReOx  Tt  D2  B21  \  Reog  T  j  ^  B22  \  Reo^  T  j  ^ 


(4.34) 


ReoJ  Tj  _  1 
Reo^  Tt  T>i 

To  obtain  measurements  in  the  x-z  plane  replace  v  by  w  and  z  replaces  y.  Also  (f>  becomes 
90  -  (/>  to  remain  consistent  with  a  right  hand  coordinate  system. 


ReoJ  TJ 


ReOf  Tt 


ReoJ  Tj 

ReOg  Tt 


(4.35) 
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J^.SA  Separation  of  Turbulence  Variables,  The  formulas  used  for  the  separation 
of  turbulent  quantities  are 


u  pu  p 


^  ^  Jl  £ 

u  pu  ^  p 


(4.36) 


where  a  =  [1  +  0.5(7  “  and  /3  =  (7  -  l)aM^.  If  p'  is  assumed  to  be  zero,  which 

may  be  problematical,  the  cross- wire  variable  may  be  separated  using  equation  4.36.  Note, 
Kistler  (7)  has  suggested  that  p'  is  proportional  to  thus  it  is  second  order.  In  addition, 
Bowersox  and  Schetz  (4)  have  experimentally  verified  the  validity  of  this  assumption  for  a 
Mach  4.0  free  mixing  layer.  The  six  equations  for  the  separated  variables  is  given  by: 


1 

0 

1 

2 

0 

0 

0 

1 

Rl 

0 

2Ro 

0 

fS^Rl 

— 2a/3 

—2aPRo 

2P'^R, 

0 

0 

Ro 

Ro 

1 

1 

0 

—a 

P  —  a 

PRo 

PRo 

0 

PRo 

—aR„ 

PRo 

PRl-a 

p 

( (#)’  1 


n 

( 

^{pvy^ 

(w) 

2 

{puY  {pvY 
pu  pu 

(puY  Tt* 

^  Ti 

I  / 

\  T,  / 


(4.37) 


These  transformations  to  Favre  averaged  N-S  equations  are  given  by  (2) 


^  =  U  =  pulp  (4.38) 

'pu"v"  —  pu'v' 


4.3.2  Single  Overheat.  If  the  hot-wire  is  operated  at  a  higher  overheat,  then 
the  sensitivity  to  the  total  temperature  fluctuations  is  minimized.  Thus  if  T/  is  small, 
then  the  SOH  technique  can  be  used  to  obtain  accurate  mass  flux  results.  The  Strong 
Reynold’s  analogy  was  developed  to  correct  for  neglecting  the  total  temperature  sensitivity 
of  the  probe.  If  the  total-temperature  fluctuations  are  assumed  negligible  (T/  =  0),  then 
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according  to  reference  (12) 


(^  +  Ro^) 

\u  u  J 


(4.39) 


Bowersox  (2)  postulated  that  when  Pr  7^  1  and  T/  0  the  temperature  fluctuations  can 
be  represented  by 

Z  =  +  (4.40) 

T  \u  u  J 

where  n  is  determined  empirically  as  ranging  from  0.0  to  1.0.  Defining  6  =  ~  1)M^, 

and  assuming  R\9  <C  1  and  <C  1  —  ^  then 


-  =  ^  [(H'  ^ 

u  1  —  0  [  'pu  ^  'pu  , 


(4.41) 


^  _  RpO  {pu)'  ^  {pv)' 
u  1  —  6  'pu 

II  ^  i3  +  (xo  r(^  ^  ^{pvy 

Ti  i-^L^  °  W' . 

with  this,  the  cross-wire  equations  can  be  reduced  to 


(4.42) 


(4.43) 


(4.44) 
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where 


Sn  =  fi^h  +  2/i5fii2mT-KuT5ii  +  g^R^rp 

Su  —  ‘^[fiBiiB2i  +  figiRmT{RuTRoBn~^ 

RvTB2l)+glKTRo] 

Sl3  “  f\B\\  +  ^fxgiRmTRvrRoB^l 
^'21  =  /I -^12  +  2/252^mr-RuT^12  +  Q^RmT 
S22  =  2[/2^5i2522  +  f292RmT{RuTRoBi2  + 

RvtB22)  +  dlRmTRo]  (4.45) 

^23  =  /I ^22  +  ‘^f292RmTRvTRoB22 
S3I  =  flf2BnBi2  +  fl92RmTRuTBii  + 
f29lRmTRuTBi2  +  9l92RmT 

S32  —  flf2{BllB22  +  B12B21)  +  fl92RmT{RuTRoBii  +  RvtB2i) 

+  f29lRmTiRuTRoBi2  +  i?i;T^22  )  +  '^9i92R^tRo 
S33  =  f\f2B2lB22  +  fl92RmTRvTRoB21  +  f29lRmT  RvtRoB22 

The  correlations  above  are  given  by: 


p  _  P-\-oi6 
^mT  —  iSf 


p  ipuyn 

Jtt  T  —  - 


ViP^y^y/W 


R 


vT 


ipvyT; 


(4.46) 


4-4  Turbulence  Transformation 

Bowersox  and  Schetz  (4)  have  shown  that  the  Reynolds  shear  stress  can  be  expressed 
in  ternas  of  the  conservative  cross-wire  variables  as 


rT  =  +Tm^(^Y 


(4.47) 


where  the  second  term  on  the  right  hand  side  has  been  shown  to  be  much  less  than  the 
first  term  for  thin  layer  type  flows.  If  the  effects  of  the  pressure  fluctuations  on  the 
cross- wire  response  are  assumed  small,  then  the  cross- wire  results  can  be  decomposed  into 
“nonconservative”  variables  (Bowersox  (3)),  i.e, 
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(p')^  (^^)^  P'^'j 

Hot  wire  anemometry  was  used  in  the  mode  of  constant  temperature.  Constant 
temperature  (resistance)  was  chosen  because  with  constant  current  anemometers  wires 
can  self  destruct  in  flows  which  have  sharp  decreases  in  velocity  (16).  The  AFIT  wind 
tunnel  has  strong  un-starting  shock  waves  which  could  possibly  burn  the  probes  out  under 
these  conditions.  In  this  study  a  constant  temperature  anemometer  with  fixed  resistance 
was  connected  to  a  1:1  bridge.  The  current  resistance  level  was  set  to  10.4  ohms.  The 
signal  of  the  hot  wire  probes  were  broken  down  into  three  components,  the  mean  value, 
the  cross-correlation,  and  the  root  mean  square  of  the  deviation. 

Mean: 


(4.48) 

V!  =  {V,-V,)/[V^{N-l)] 

(4.49) 

VlVi  =  {v;  *  V^)/[{N  -  1)  *  *  F2] 

(4.50) 

Where  N  is  the  number  of  sampled  points,  and  i  is  the  wire  number. 
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V,  Results  and  Discussion 


A  characterization  of  the  boundary  layer  and  its  turbulent  structure  was  obtained 
using  the  Pitot  probe,  cone  static  probe,  and  the  cross  wires  at  a  location  44  cm  from 
the  nozzle  throat  to  document  the  flow  before  the  expansion  region.  In  this  way,  when 
taken  in  concert  with  the  same  measurements  (Pitot,  cone  static,  and  cross  wires)  made  in 
the  favorable  pressure  gradient  region  (71.5  cm  downstream  from  the  nozzle  throat),  the 
effects  of  the  changes  in  the  pressure  on  the  turbulent  flow  structure  can  be  ascertained. 


5,1  Shadowgraph  and  Schlieren  Photography 


Figure  5.2  Test  Section  Shadowgraph  Image 
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Figure  5.3  Test  Section  Mach  Profile  x  =  71.5  cm 

In  Figure  5.1  is  presented  a  schlieren  image  in  the  vicinity  of  the  expansion  fan.  Note 
that  the  boundary  layer  has  not  separated,  the  thickness  of  the  boundary  layer  is  increasing 
in  the  streamwise  direction,  and  the  presence  of  turbulent  eddies.  The  shadowgraph  image 
in  Figure  5.2  is  presented,  but  is  of  little  value  in  an  expansion  region  since  the  mean  density 
gradients  are  rather  small.  In  Figure  5.3  is  presented  the  profile  of  the  Mach  number  across 
the  wind  tunnel  as  determined  using  Pitot  and  cone-static  probes.  Figure  5.4  is  a  plot  of 
the  PniPtoo  profile  at  x  =  71.5  cm.  The  jumps  in  pressure  are  located  in  the  vicinity  of 
the  expansion  fans. 

5.2  Conventional  Data 

Detailed  Pitot  and  cone-static  measurements  are  obtained  along  the  tunnel  centerline 
at  X  =  44  cm  and  71.5  cm.  Figure  5.5  shows  a  rough  outline  of  the  measurement  locations, 
where  1  is  the  x  =  44  cm  location  and  2  is  the  x  =  71.5  cm  location.  The  latter  station 
is  located  in  the  vicinity  of  the  expansion.  Variation  of  the  mean-flow  profiles  in  the 
streamwise  direction  highlights  the  effect  of  the  pressure  gradient;  the  mean  flow  accelerates 
under  the  influence  of  a  favorable  pressure  gradient  (Figure  5.3).  The  boundary  layer  edge, 
at  the  44  cm  location,  is  defined  as  the  point  at  which  the  mean  velocity  reached  99%  of 
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Figure  5.4  Test  Section  Pressure  Profile  PnIPtoo  x  ==  71.5  cm 

the  mean  value  at  the  tunnel  centerline.  Because  the  flow  in  the  expansion  section  is  not 
fully  developed,  the  99%  criteria  is  based  on  the  highest  Mach  number  peak  near  the  wall. 
This  peak  Mach  number  is  3.0.  The  boundary  layer  thickness  at  x  =  44  cm  is  estimated 
at  0.83  cm  and  1.17  cm  downstream  (x  =  71.5  cm).  The  freestream  Mach  number  is  2.9 
upstream.  Measurements  acquired  as  part  of  the  current  study  extend  to  within  0.06 

of  the  wall  upstream  and  yjS  ^  0.09  downstream. 


Figure  5.5  Test  Section 


5,3  Repeatability  of  Data 

One  difficulty  with  experimental  work  is  that  the  data  must  be  measured  in  a  stable 
environment  where  the  variable  of  the  conditions  do  not  change.  The  precautions  men¬ 
tioned  earlier  in  Chapter  III  highlight  some  of  the  steps  taken  to  insure  the  equipment  is 
stable.  As  a  means  of  investigating  the  repeatability  of  the  data,  several  duplicate  runs 
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were  conducted.  In  particular,  the  profiles  of  the  flow  in  Figure  5.6  were  taken  on  two 
different  days  with  the  same  hot-wire  probe.  The  profiles  of  the  flow  in  Figure  5.7  are  from 
the  second  channel  of  the  UW  x-wire  probe.  Through  comparison  of  flow  profiles  obtained 
using  alternative  means  (i.e.  conventional  and  hot  wire  methods)  a  comparison  was  made 
of  Mach  profile  shape  and  voltage  cross-correlation  from  the  cross- wires.  Figure  5.8  shows 
the  Mach  profile  through  the  boundary  layer  and  Figure  5.9  shows  the  downstream  shear 
profile  from  the  SOH  data.  Again,  the  close  correlation  of  the  profile  shapes  demonstrates 
the  repeatability  of  the  data.  To  confirm  the  consistency  of  the  hot  wire  data,  the  mean 


Figure  5.6  UW  RMS  Voltage  Wire  l(x  =  44  cm) 

flow  data  profile  shapes  are  compared.  The  mean  flow  voltages  are  plotted  in  Figure  5.10. 
The  close  agreement  in  shape  of  the  five  plots  confirms  the  repeatability  of  the  runs. 

5.4  Hot  Wire  Constants 

The  constants  in  the  equations  for  reduction  of  the  SOH  data  are  derived  from 
the  MOH  data.  The  terms  Rvt  ^.nd  R^jt  were  determined  by  averaging  the  values 
calculated  from  different  MOH  runs.  They  were  determined  to  be  0.6,  -0.5,  and  0.15.  The 
variable  k,  was  determined  by  matching  the  SOH  and  MOH  data  by  adjusting  k  until  the 
data  aligned.  This  value  was  highly  dependent  on  the  SOH  ratio  and  can  be  shown  that 
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Figure  5.7  UW  RMS  Voltage  Wire  2(x  =  44  cm) 


Mach  Number 


Figure  5.8  Mach  Profile  (x  =  44  cm) 


K  =  f{OH R).  The  overheat  ratios  (OHR)  force  a  sensitivity  to  flow  properties  that  respond 
to  factors  in  the  Strong  Reynolds  analogy  equation.  The  Strong  Reynolds  analogy  assumes 
K  ^  1(11)-  This  applies  to  flows  where  the  sensitivities,  f  and  g,  are  not  coincidental  and 
are  widely  separated,  i.e.,  overheat  levels  below  approximately  1.8.  The  sensitivities  are 
mentioned  in  Chapter  4  and  are  related  to  the  system  fluctuations  by 


V  'pu 


+  9 


V 

% 


(5.1) 


An  approximate  graph  of  the  behavior  of  f  and  g  looks  like  Figure  5.11  and  ^  is  usually 


Figure  5.11  Hot  Wire  Probe  Sensitivities 


assumed  small.  To  solve  for  k  experimentally  multiple  overheat  values  are  used  and  flow 
properties  are  derived  as  in  chapter  4.  In  this  way  the  data  is  made  independent  of  overheat 
level.  For  overheat  values  larger  than  1.9,  k  has  a  more  gradual  influence. 

For  the  SOH  UV  data  the  overheat  level  was  approximately  1.6,  and  a  value  of 
K  Ri  1.0  is  appropriate,  as  postulated  in  Morkovin’s  work.  This  aligned  the  UV  data  with 
the  MOH  data.  Figure  5.12  represents  k  =  0,  Figure  5.13  represents  k  =  0.4,  and  Figure 
5.14  represents  k  =  1.0. 

The  UW  SOH  data  on  the  other  hand  was  found  to  align  best  at  k  =  0.  The  UW  data 
overheat  ratio  was  approximately  2.2,  causing  the  data  to  be  less  sensitive  to  k  changes. 
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Figure  5.15  shows  UW  data  with  k  =  1.0  and  Figure  5.16  shows  when  k  =  0.0.  The 
K  =  0.0  data  aligns  best  with  the  multiple  overheat  data  and  removes  the  disparity  from 
the  U  component  in  the  UV  and  UW  data.  When  k  was  set  to  either  1.0  or  0.4  for 


Figure  5.13  UV  Probe  Turbulence  Intensities  (k  =  .4  x  =  44  cm) 

both  channels  the  U-component  differed  by  0.02.  The  sensitivity  disparity  between  f  and  g 
are  accounted  for  by  k  =  1.  By  increasing  the  overheat  ratios  to  be  near  2.0  the  disparity 


5-8 


Figure  5.15  UW  Probe  Turbulence  Intensities  ( 


between  f  and  g  sensitivities  can  be  removed  when  k  is  set  to  zero,  i.e.  less  temperature 
sensitivity  in  equation  4.40  (T’/Tt  =  0). 

5.5  Zero  Pressure  Gradient  Region 

The  turbulence  intensities  represented  in  Figure  5.14  and  Figure  5.16  of  the  UV  and 
UW  wires,  respectively  are  the  quantities  measured  directly  by  the  hot-wires.  The  relative 
magnitude  of  this  term  is  the  largest  of  the  intensity  terms  in  the  boundary  layer.  This 
term  also  appears  to  be  asymptotically  approaching  a  much  lower  value  as  it  nears  the 
wall.  The  terms  ((pu)7(pu)),„,  and  {{pw)' j {pu))rm3  both  appear  to  be  approaching  the 
wall  at  a  linear  rate.  The  multiple  overheat  data  shows  good  agreement  with  the  single 
overheat  data. 

The  separated  turbulence  intensities  follow  a  similar  pattern.  The  Figures  5.17  and 
5.18  illustrate  the  separated  turbulence  intensities.  The  components  of  u' ,  v',  w'  all  follow 
a  nearly  linear  slope  through  the  boundary  layer.  Near  the  wall  all  three  terms  collapse 
within  the  same  region  of  both  graphs.  Additionally,  when  the  terms  reach  the  edge  of  the 
boundary  layer  {v'/u\ms,  {w'lu\,ns  and  {p' j p)rms  collapse  into  the  same  region,  whereas 
the  term  {u'/u)rm3  is  near  zero  in  the  freestream.  The  fluctuating  term  {p'/p)rms  is  the 
largest  fluctuating  component  through  most  of  the  boundary  layer. 

The  component  of  shear  in  the  x-y  plane  is  predominantly  composed  of  the  up'v' I'piP 
term.  Figure  5.19.  Within  the  majority  of  the  boundary  layer,  neglecting  the  boundary 
layer  region  below  yjS  k,  0.2  and  the  top  10%  of  the  boundary  layer  edge,  up'v'j'piP 
is  approximately  67%  of  the  total  shear.  Neglecting  this  term  would  be  neglecting  the 
majority  of  the  shear.  The  terms  v'  and  p'  are  the  larger  fluctuating  components  in  the 
boundary  layer  and  contribute  to  create  the  larger  term  up'v'  j'piP.  Similarly,  the  terms 
w'  and  p'  contribute  to  the  up'w'j'plP  term  being  the  largest  contribution  to  the  term. 
Figure  5.20.  Note  that  the  maximum  value  of  T^y  is  twice  as  large  the  maximum  of 
Furthermore,  the  term  is  only  significant  in  the  lower  half  of  the  boundary  layer. 
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Figure  5.18  UW  Probe  Separated  Turbulence  Intensities  (x  =  44  cm) 
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Figure  5.19  UV  Probe  Shear  (x  =  44  cm) 
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The  UV  Reynolds  shear  derived  from  the  MOH  data  is  in  close  agreement  with  the 
measured  shear  stress.  The  MOH  derived  values  of  the  Favre  shear 

rjy  \Favri=  pu'v"  ^  (5.2) 

estimates  the  first  term  of  the  Reynolds  shear  well  also  showing  the  close  agreement  between 
the  SOH  and  MOH  data.  Morkovin  postulated  that  the  velocity  density  correlation, 

up'v' 

-  «  1 
pu'v' 

where  we  have  found  this  quantity  to  be  approximately  2  though  60  to  70%  of  the  boundary 
layer. 

5.6  Favorable  Pressure  Gradient  Region 

The  raw  turbulence  intensities  {{pu)' /{pu))rms,  etc.  in  the  expansion  region  are  still 
dominated  by  the  term  ((pu)' l{pu))rms-,  Figure  5.21  and  Figure  5.22.  The  freestream 
value  of  the  UW  probe  is  0.02  whereas  the  UV  probe  has  a  value  of  approximately  0.03  in 
the  freestream.  The  components  ((pw)' / (pu))rms  and  ((pv)' j (pu))rms  both  follow  a  nearly 
linear  progression  from  the  freestream  value  down  to  maximum  values  of  0.06  and  0.1 
respectively.  The  component  ((pu)' j (pu))rms  follows  a  ’’fuller”  curved  path  which  curves 
back  towards  the  value  of  the  other  two  terms  as  it  approaches  the  lowest  measured  station. 

The  separated  turbulence  intensities  provide  an  insight  as  to  which  terms  are  caus¬ 
ing  this  characteristic,  Figure  5.23  and  Figure  5.24.  Again  the  predominant  fluctuation 
quantity  is  the  density,  (p' j p)rms-  As  the  flow  approaches  the  wall,  below  yjb  k,  0.2,  the 
components  of  u'  and  w'  become  the  predominant  terms,  u',  v' ,  and  w'  vary  linearly  across 
the  boundary  layer,  while  the  term  (p'  j p)rms  has  the  characteristic  shape  of  the  turbulence 
intensity  term  ((pu)' j (pu))rms- 

From  these  observations  and  experience  with  the  upstream  (44  cm)  data  one  could 
conclude  that  the  shear  stress  would  predominantly  be  driven  by  the  density  fluctuation 
component.  The  component  of  the  shear.  Figure  5.25  and  Figure  5.26,  up'v'lplf  is 
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Figure  5.20  UW  Probe  Shear  (x  ==  44  cm) 
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Figure  5.21  UV  Probe  Turbulence  Intensities  (x  =  71.5  cm) 
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the  dominant  term,  exhibiting  the  same  behavior,  with  67%  of  the  shear  residing  in  this 
term  throughout  most  of  the  boundary  layer.  The  shear  is  less  pronounced,  but  the 
density  term  up'w'  1'^^  accounts  for  at  least  50%  of  the  shear  in  this  plane.  The  MOH 
Reynolds  shear  estimation  for  the  x-y  plane  tracked  the  SOH  term  up'v' IJTU?  throughout 
the  boundary  layer.  This  term  accounts  for  67%  of  the  total  shear.  The  Favre  on  the 
other  hand  continued  to  track  the  pu'v”  jpu^  term.  Given  the  predominance  of  the  mean 
flow  in  the  streamwise  direction,  shear  stress  terms  such  as  -vp'u'l^'^  and  -wp'u'lpu'^, 
are  insignificant.  Similar  to  the  upstream  flow  the  velocity-density  correlation  relation  of 
Morkovin  is  also  found  to  be  ~  2.0. 

Comparison  of  the  turbulence  intensities  from  the  44  cm  upstream  test  station  to 
the  71.5  cm  test  station  shows  subtle  changes  in  the  character  of  the  flow  (  Figures  5.14, 
5.21,  5.16  and  5.22).  The  component  {{pv)' ! {pu)\mz  increased  slightly  in  the  downstream 
location.  The  freestream  value  of  both  of  the  u  and  v  components  increased  from  ap¬ 
proximately  0.02  upstream  to  0.03  downstream.  This  change  accounts  for  the  increase  of 
the  {{pv)' ! {pu)\ms  data.  The  UW  wire  shows  a  similar  trend.  A  slight  increase  in  the 
{{pw)' f  {pu))rms  component  can  be  seen  in  the  freestream.  Additionally,  the  peak  value  of 
((pu)' I {pu))rms  has  decreased  downstream,  but  it  also  has  a  ’’fuller”  profile.  The  separated 
turbulent  intensities  show  similar  characteristics.  The  maximum  value  of  v'  increases  from 
a  peak  value  upstream  of  «  0.085  to  «  0.1  downstream.  The  components  of  turbulence  for 
u'  remains  relatively  constant,  with  a  very  slight  decrease  in  the  peak  value  downstream. 
The  slope  of  the  u'  and  v'  component  through  the  boundary  layer  does  not  change.  The 
peak  of  the  density  fluctuation  decreases. 

At  the  top  and  bottom  10%  of  the  boundary  layer  (Figures  5.19,  5.25,  5.20  and 
5.26),  the  proportionality  of  the  shear  components  is  approximately  50/50.  Throughout 
the  remainder  of  the  boundary  layer  the  proportionality  is  split  with  33%  in  the  first  term 
and  67%  in  the  second  term.  It  is  curious  to  see  that  the  first  contribution  to  the  Reynolds 
stress  appears  to  be  increasing  asymptotically  as  it  approaches  the  lowest  measured  station. 
The  second  term  seems  to  die  off  rapidly  after  passing  the  point  oi  yjb  ~  0.2.  This  last 
observation  is  not  as  apparent  in  the  upstream  graph. 
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The  reduction  of  the  shear  into  the  separate  Reynolds  stress  components  assumes 
that  the  second  term  in  the  conservative  cross- wire  shear  equation  is  small. 


P  P 


(5.3) 


This  can  be  verified  by  looking  at  the  graph  of  the  shear  contributions,  Figure  5.27  and 
Figure  5.28.  The  second  term  is  essentially  a  third  order  term  because  v  is  on  the  order  of 
u',  causing  this  term  to  be  small.  So  although  the  pressure  fluctuations  must  be  assumed 
small  to  separate  out  the  terms  in  the  Reynolds  stresses  equation,  i.e. 


tJ-  =  -pu\u'j  -  Uip'u'j  -  Ujp'u'f  -  p'u'iU'j  (5.4) 

the  value  of  the  Reynolds  stress  can  be  assumed  to  be  accurate  regardless  of  the  pressure 
fluctuation  assumption. 


5. 7  Heat  Flux 

Few  studies  have  documented  the  turbulent  heat  flux  behavior.  The  use  of  MOH 
allows  the  separation  of  the  heat  flux  terms  necessary  for  solution  of  the  energy  equation. 
Figure  5.29  shows  the  heat  flux  terms  from  the  UV  cross  wire  and  Figure  5.30  shows  the 
values  obtained  from  the  UW  cross  wire. 
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VI.  Conclusions  and  Recommendations 


6.1  Conclusions 

The  importance  of  compressibility  and  the  effect  of  pressure  gradient  on  the  turbu¬ 
lent  flow  structure  in  a  Mach  2.9  boundary  layer  were  investigated  experimentally.  Con¬ 
sistently,  terms  involving  density  fluctuations  are  shown  to  be  important  in  both  the  zero 
and  favorable  pressure  gradient  regions. 

This  effect  manifests  itself  in  the  UiP'u'j  term  appearing  in  the  RANS,  accounting 
for  as  much  as  67%  of  the  total  shear  stress  (irrespective  of  pressure  gradient).  The  first 
two  terms  in  do  not  behave  in  a  linear  manner  but  combine  to  create  an  overall  linear 
total  through  the  boundary  layer.  The  shear,  is  remarkably  linear  through  most  of 
the  boundary  layer  with  a  slope  upstream  of  approximately  —182  and  downstream  of 
—  109.  The  upstream  value  is  approximately  zero  in  the  freestream  and  estimated  as  .0008 
downstream. 

In  addition,  the  fluctuating  components  in  the  energy  equation  may  be  of  similar 
importance.  The  terms  in  the  energy  equation  hoP'u'i  and  UiP'h'g  may  be  significant  because 
they  contain  the  density  fluctuation.  The  FANS  equations  again  do  not  contain  terms 
involving  explicit  mention  of  the  density  fluctuation  component  in  the  energy  equation. 

The  density  fluctuation  terms  in  the  stress  calculations  are  a  significant  component 
of  the  shear  stress.  The  lack  of  an  explicit  term  in  the  FANS  equations  to  account  for 
this  fluctuation  demonstrates  the  difficulty  of  modeling  density  fluctuations  in  the  FANS 
based  CFD  solutions.  CFD  methods  should  account  for  compressibility  effects  explicitly  in 
order  to  more  accurately  model  the  high  speed  turbulent  flow  regime.  Through  the  FANS 
form  of  the  N-S  equation  is  more  straightforward  to  model  numerically,  most  experimental 
techniques  measure  RANS-type  terms.  Thus,  in  flow  situations  where  the  importance  of 
compressibility  and  pressure  gradient  can  be  identified,  the  incorporation  of  this  informa¬ 
tion  into  turbulent  models  is  more  suited  for  RANS  applications.  Thus  if  the  interest  is  to 
develop  more  accurate  closure  models  based  on  this  data,  a  RANS  approach  may  be  most 
useful. 
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The  assumption  that  k  =  1.0  suggested  by  in  the  Strong  Reynolds  analogy  is  not 
universally  applicable.  Overheat  ratios  of  about  1.6  are  likely  candidates  for  the  use  of 
K  =  1.0.  On  the  other  hand  this  investigation  has  found  if  the  overheat  ratio  is  near  2.0  a 
K  0.0  is  more  appropriate.  The  study  overheat  ratios  of  1.6  for  the  UV  data  matched  the 
MOH  data  best  when  the  value  of  1.0  was  used.  However,  when  the  overheat  ratios  were 

2.2  the  value  of  k  of  0.0  matched  the  MOH  best.  The  Strong  Reynolds  analogy  is  correct 
that  if  a  moderate  overheat  ratio  is  used  that  1.0  is  an  ideal  number  for  matching  SOH 
and  MOH  data,  but  as  the  overheat  level  becomes  higher  the  sensitivity  of  the  data  to  k 
decreases,  therefore  use  of  the  Strong  Reynolds  analogy  was  only  necessary  at  the  lower 
overheat  ratios. 

6.2  Recommendations 

To  further  quantify  the  effect  of  a  favorable  pressure  gradient,  measurements  should 
be  obtained  of  a  zero  pressure  gradient  boundary  layer  at  the  same  location  as  where  the 
favorable  pressure  gradient  measurements  were  conducted  (71.5  cm).  In  this  way,  changes 
in  turbulence  intensities  and  shear  stress  due  only  to  the  evolution  of  the  boundary  layer 
in  the  absence  of  pressure  gradients  could  be  extracted,  and  a  truer  measure  of  the  role  of 
the  pressure  gradient  obtained. 

LDV  (Laser  Doppler  Velocimetry)  measurements  of  same  flow  should  be  conducted 
in  order  to  assess  the  effect  of  the  hot-wire  probes  on  the  flow  field. 

The  use  of  higher  overheat  ratios  in  the  SOH  data  portion  of  this  experiment  would 
remove  the  need  for  a  k  other  than  zero.  The  data  of  the  SOH  would  match  the  MOH  data 
more  closely  without  data  matching  between  the  MOH  and  the  SOH  data.  Because  the 
lower  overheat  ratios  provided  a  better  picture  of  the  temperature  fluctuations,  overheat 
values  and  the  range  of  the  MOH  should  have  spanned  the  overheat  levels  more  fully,  i.e. 
from  1.1  up  to  2.2  if  possible. 


6-2 


Appendix  A.  Conventional  Probe  Data  Plots 


Figure  A.l  Upstream  Pressure  Profile  Pn! Pt\  x  =  44  cm 


Figure  A. 2  Test  Section  Pressure  Profile  Pt2l Pn  x  =  71.5  cm 


Appendix  C.  SOH  Cross  Correlations 
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Figure  C.l  UV  Cross  Correlation  Voltage(x  =  44cm) 
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Figure  C.2  UW  Cross  Correlation  Voltage(x  =  44cm) 


/oltage  Wire  l(x  =  71.5  cm) 


Appendix  E.  Turbulence  Intensities  and  Separated  Turbulence  Intensities 
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Figure  E.l  UV  Probe  Turbulence  Intensities  (x  =  44  cm) 
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Figure  E.2  UV  Probe  Turbulence  Intensities  (x  =  71.5  cm) 
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Figure  E.4  UW  Probe  Turbulence  Intensities  (x 


Appendix  F.  Shear  Stress 
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Figure  F.l  UV  Probe  Shear  (x  =  44  cm) 
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Figure  F.2  UV  Probe  Shear  (x  =  71.5  cm) 


Figure  F.4  UW  Probe  Shear  (x 


F.2  FANS 


Figure  F.5  UV  Probe  Shear  (x  =  44  cm) 
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Figure  F.6  UV  Probe  Shear  (x  =  71.5  cm) 


Figure  F.12  UW  Probe  Shear  (x  =  71.5  cm) 


Appendix  G.  Compressible  Shear  Stress 
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Figure  G.l  UV  Probe  Shear  (x  =  44  cm) 


Figure  G.2  UV  Probe  Shear  (x  =  71.5  cm) 


Appendix  H,  Separated  Favre  Turbulence  Intensities 
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Figure  H.l  UV  Separated  Turbulence  Intensities  (x  =  44  cm) 
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Figure  H.2  UV  Separated  Turbulence  Intensities  (x  =  71.5  cm) 
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Figure  1.2  UW  Mass  Flux- Total  Temperature  Correlations  (x  =  44  cm) 


Appendix  J.  Velocity-Density  and  Velocity-Velocity  Correlations 
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Figure  J.3  UV  Probe  Turbulence  Correlations  (x  =  44  cm)  Velocity-Density  Correlations 


Figure  J.4  UW  Probe  Turbulence  Correlations  (x  =  44  cm)  Velocity-Density  Correlations 
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Figure  J.5  UV  Probe  Turbulence  Correlations  (x  ==  71.5  cm)  Velocity- Density 
Correlations 
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Figure  J.6  UW  Probe  Turbulence  Correlations  (x  =  71.5  cm)  Velocity-Density 
Correlations 
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Appendix  K,  Procedure 


K.l  Mean  Flow  Data 

The  initial  reduction  of  the  data  was  taking  the  raw  Multipro  file,  which  was  binary 
and  processing  it  into  an  ASCII  file  using  a  supplied  program  by  the  manufacturer.  Al¬ 
though  more  usable  it  was  in  a  convoluted  form.  The  essential  data  contained  in  these 
files  is  the  number  of  data  points,  the  Multipro  internal  offset  value,  the  Multipro  internal 
slope  value,  the  user  entered  offset  and  the  user  entered  slope.  The  internal  file  vlaues  for 
the  user  offset  and  slope  was  not  used  and  the  calibration  values  were  used  in  the  data 
reduction.  Once  these  values  had  been  stripped  from  the  file,  the  data  is  all  that  remained. 
A  code  was  written  that  took  the  original  ASCII  file  and  constructed  a  more  compact  and 
usable  one  only  containing  the  data  and  the  above  mentioned  values.  Next,  the  data  up¬ 
stream  was  taken  at  1.27  cm  increments  to  allow  high  resolution  of  the  flow.  This  required 
five  passes  through  the  flow  to  measure  the  complete  cross  section.  The  data  had  to  be 
spliced  together  at  the  junctures  of  commonality.  Within  the  aforementioned  code  was  a 
filter  to  detect  when  the  LVDT  had  started  and  had  stopped  traversing  the  flow,  at  these 
points  the  data  was  spliced  together. 

K,2  SOH  Data 

The  single  overheat  data  began  essentially  the  same  as  the  mean  flow  data  reduction. 
The  conversion  of  the  Multipro  files  to  a  usable  ASCII  file  was  done  by  using  the  code 
AsctoAscii.c.  The  data  for  the  SOH  was  taken  in  one  traverse  step  to  a  distance  of 
approximately  two  boundary  layer  thicknesses.  This  distance  was  estimated  from  the 
shadowgraphs  and  schlieren  photographs.  Due  to  traverse  equipment  restrictions  and  the 
time  limit  of  21  seconds  for  the  run  time  of  the  tunnel  the  travel  distance  was  usually 
restricted  to  below  2.54  cm.  The  calibration  data  was  obtained,  generally  on  the  same 
run  as  the  data  by  varying  the  pressure  in  the  wind  tunnel  after  the  probe  had  stopped 
traversing.  In  this  manner  data  acquired  on  a  run  was  guaranteed  to  have  calibration  data 
available.  A  script  file  was  written  on  the  Sparc  20’s  to  execute  the  various  data  reduction 
codes.  First  was  the  ASCII  conversion,  second  the  LVDT  data  was  converted  into  inches 
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using  the  Convert.c  code.  After  the  LVDT  data  was  in  inches  the  code  LVDTcountandcut.c 
was  used  to  capture  the  intervals  where  the  LVDT  was  moving  and  count  there  position 
in  the  file  and  also  cut  this  data  out  of  the  overall  LVDT  data.  The  count,  positions, 
obtained  in  this  manner  were  used  to  cut  the  other  files  containing  the  pressure,  and  two 
channel  of  hot  wire  data.  The  code  used  for  this  was  the  Cutfile.c  code.  The  Cutfile.c  code 
had  an  additional  function  of  cutting  the  calibration  data  off  of  the  end  of  the  data  files. 
After  the  files  were  cut  they  were  converted  to  voltages  or  pressures  using  the  Convert.c 
program.  By  cutting  the  three  files,  pressure  ,  channell  hot  wire,  channel2  hot  wire,  first, 
before  Convert.c,  the  reduction  time  of  the  data  was  cut  in  half.  The  calibration  files  were 
then  converted  to  the  appropriate  units.  Once  the  calibration  files  were  converted  they 
were  converted  into  files  intended  for  MSHeaR  (2). 

The  format  for  the  MSHeaR  calibration  file  is: 

Voltage  Pressure(P(i)  Temperature(Tt)  Mach  Pn! Pn 

After  this  point  the  calibration  files  were  edited  to  remove  the  tunnel  unstart  data  from 
the  files,  this  is  characterized  by  large  fluctuations  in  data  values.  The  calibration  data 
was  very  near  linear,  as  it  should  be  in  this  Mach  range.  Next  the  construction  of  a 
mean  flow  data  file  that  matched  the  data  point  locations  of  the  SOH  data.  The  program 
Interpolate. c  was  used  to  linearly  interpolate  the  original  mean  flow  data  so  that  the  y 
coordinates  matched  the  SOH  data.  This  was  another  file  required  by  MSHeaR.  The  data 
from  the  two  channels  of  hot  wire  data  were  cross  correlated  and  there  mean  values  and 
root  mean  square  values  were  found  using  the  code  Xwk.c  written  by  Mike  Senseney  and 
Raymond  Miller.  The  original  routine  was  written  for  Matlab  by  Mike  Senseney  and  was 
subsequently  rewritten  in  C  by  Raymond  Miller.  This  resulted  in  a  conservative  estimate 
of  a  three  time  speed  increase  in  the  code  execution.  Next  the  code  MSHeaR  was  run  to 
obtain  the  initial  turbulence  results.  Last  the  code  written  by  Jon  Dotter  and  modified 
by  Raymond  Miller  was  used  to  separate  turbulent  variable  terms.  The  code  was  called 

separt.f 
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K.3  MOH  Data 


Data  for  the  multiple  overheat  was  much  harder  to  reduce  than  either  of  the  two  other 
data  types.  A  script  file  was  also  written  to  execute  most  of  these  steps  automatically.  First 
you  have  a  separate  calibration  data  set  that  may  or  may  not  have  the  resistance  levels  in 
the  same  order  as  your  data  file  and  second,  filtering  the  data  is  very  time  consuming.  The 
data  is  initially  broken  down  the  same  as  the  other  data,  conversion  to  a  usable  ASCII  form 
using  AsctoAscii.c.  The  data  is  then  converted  to  its  native  form,  voltage,  etc.  because  the 
data  had  to  be  viewed  to  set  various  data  cuting  criteria.  The  method  in  which  the  data 
was  triggered  made  it  no  longer  crucial  to  base  data  reduction  on  the  LVDT,  but  on  one  of 
the  cross  wire  channels.  The  two  channels  of  data  from  the  hot  wires  were  viewed  and  the 
more  stable,  least  wildly  oscillatory,  was  used  for  the  cut  indexing.  This  channel  was  cut 
on  two  criteria,  first  if  there  was  a  false  jump  then  a  skip  first  jump  switch  was  set,  second 
the  smallest  jump  between  overheat  steps  was  set  as  teh  second  criteria.  These  jumps 
ranged  in  value  from  0.6  volts  up  to  4.0  volts  typically.  Usually  the  larger  was  used  and 
then  cut  down  until  a  satisfactory  cut  was  made.  Another  adjustment  that  could  be  made 
within  Channelcut.c  was  the  number  of  data  points  to  not  use  at  the  end  of  the  step  of 
data.  For  example  if  the  jump  size  was  set  to  4.0  and  the  data  had  large  jumps  in  it  at  the 
end  of  the  step  then  the  cut  at  the  end  of  the  data  would  often  ignore  350  points  and  grab 
the  2056  points  that  proceed  this  area.  In  this  way  adjustments  could  be  made  as  to  where 
on  the  data  step,  overheat  ratio,  the  data  was  taken  from.  Typically  the  lower  the  criteria 
the  earlier  the  initial  cut  location.  By  trial  and  error  the  calibration  files  and  the  data  file 
that  was  used  as  a  key  had  to  be  cut  ’’cleanly”.  Cleanly  meaning  no  wild  oscillations  of 
transition  inherent  in  the  data.  Once  the  two  channels  used  as  keys  were  satisfactory  the 
other  channels  had  to  be  cut  using  Cutfile.c,  same  name  different  algorithm,  to  match  the 
key  channel.  These  files  had  to  also  be  viewed  to  see  if  further  adjustments  needed  to  be 
made  to  the  cut  criteria.  The  LVDT  data  was  sampled  at  a  lower  rate  than  the  pressure 
or  hot  wires  and  therefore  had  to  be  expanded  to  meet  there  number  of  data  points,  this 
was  accomplished  using  the  code  Expand. c  which  takes  the  frequencies  of  each  channel  and 
expands  the  file  to  match.  After  the  LVDT  data  is  expanded  it  is  cut  just  like  the  other 
files. 
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Once  all  files  have  been  cut  then  the  calibration  files  must  be  sorted.  Calibration  files 
were  captured  on  separate  runs  from  the  data  runs.  Sometimes  the  sequence  of  resistors 
was  in  a  different  order  from  the  data  sequence.  The  code  to  separate  the  calibration  data 
allowed  such  reordering.  The  code  Calibrate. c  constructed  the  calibration  files  needed  by 
MSHeaR. 

The  calibration  files  are  now  ready  and  the  interpolated  mean  flow  data  is  needed. 
This  code  is  the  same  as  above  and  does  the  same  linear  interpolation  of  data  to  match 
the  y  locations.  The  code  MSHeaR  is  run  at  this  point  and  the  function,  either  option  7 
or  8  is  used  to  evaluated  the  MOH  data,  8  was  used  for  all  data  in  this  thesis. 

A  bug  when  running  MSHeaR  on  Unix  systems  is  that  the  MOH  separation  routine 
needs  a  space  between  the  title  of  the  output  file  of  option  8  inorder  to  run  option  11  on 
all  data  points.  Option  11  provided  MOH  separtion  of  variables. 

K.4  Midway 

The  results  in  this  section  are  not  valid  because  the  mean  flow  data  was  not  available. 
If  an  adjustment  to  the  mean  data  or  measurements  are  taken  later  this  data  may  be  highly 
useful  in  tracking  the  flow  properties  in  the  wind  tunnel.  It  has  been  suggested  that  a 
delta  shift  of  the  44  cm  location  mean  flow  data  may  render  this  data  useful.  Intermediate 
turbulence  measurements  are  taken  53.5  cm  from  the  nozzle,  between  the  two  original 
upstream  station,  44  cm,  and  the  downstream  station,  71.5  cm.  The  measurements  are 
obtained  in  a  zero  pressure  gradient  flat  plate  section.  The  only  terms  measured  are  the 
UW  single  overheat  readings.  They  provide  some  insight  into  the  nature  and  progression 
of  the  flow.  In  order  to  provide  direct  comparison  with  the  44  cm  location,  data  is  non- 
dimensionalized  using  the  the  conventional  probe  values  from  the  44  cm  location,  i.e.  Mach 
number,  etc.  The  turbulence  intensities,  separated  turbulence  intensities  and  the  shear  all 
increase  dramatically.  Figures  K.l  K.2  and  K.3.  The  component  of  the  shear  in  the  xz 
plane  increased  dramatically.  The  shear,  has  increased  from  a  maximum  of  0.003  to 
a  maximum  of  0.009,  a  3  times  increase  in  magnitude!  The  UW  data  at  the  53.5  cm 
cite  resembles  the  UV  data  at  the  44  cm  location  in  magnitude  and  character.  Still  it  is 
quit  clear  that  the  density  fluctuation  term  is  the  largest  term  present  according  to  the 
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Table  K.l  Favre  terms  from  multiple  overheat  data 


y 

Run  009 

Run  010 

0.840E-01 

0.221E-02 

-0.102E+00 

0.301E+00 

0.145E-02 

0.289E-03 

0.517E+00 

0.678E-03 

-0.173E-03 

0.734E+00 

0.499E-03 

-0.361E-02 

0.950E+00 

0.129E-03 

-0.137E-02 

Table  K.2  Favre  terms  from  multiple  overheat  data 


y 

on 

012 

0.119E+00 

0.209E-02 

- 

0.424E+00 

0.360E-03 

0.155E-02 

0.728E+00 

0.155E-04 

0.439E-04 

0.103E+01 

-0.297E-04 

-0.836E-06 

0.134E+01 

-0.200E-04 

- 

previous  description,  approximately  67%  of  the  magnitude  of  the  total  shear  is  accounted 
for  in  this  term.  Comparison  of  the  turbulence  intensities  from  the  53.5  cm  location  to 
the  71.5  cm  location  shows  more  dramatic  changes.  The  53.5  cm  data  suggest  that  the 
boundary  layer  is  still  expanding  and  the  turbulence  is  getting  stronger.  The  expansion 
section  appears  to  damp  out  the  turbulence  and  smooth  the  flow.  The  turbulence  intensity 
{{pu)' ! {pu))rms  has  decreased  from  a  peak  value  of  0.14  to  0.10  when  traveling  through 
the  expansion  ramp.  Recall  that  the  53.5  cm  location  is  non-dimensionlized  with  the 
44  cm  data  freestream  values.  The  separated  values  show  a  similar  trend.  An  increase 
in  magnitude  from  the  44  cm  section,  but  then  a  sharp  decrease  in  the  71.5  cm  section 
when  exposed  to  the  expansion  ramp.  Once  again,  this  demonstrates  a  favorable  pressure 
gradient  damps  the  turbulence  in  a  flow. 

The  comparison  of  the  shear  for  the  53.5  cm  section  to  the  71.5  cm  section  shows  a 
more  drastic  decrease  in  the  shear  value  from  a  maximum  of  almost  0.009  to  approximately 
0.0013,  almost  a  7  fold  decrease  in  spanwise  turbulence. 


K.5  Tables  of  Favre  and  Reynolds  Shear  Terms 
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Figure  K.3  UW  Probe  Shear  (x  =  53.5  cm) 


Table  K.3  Reynolds  terms  from  multiple  overheat  data 


y 

on 

012 

0.119E+00 

0.373E-02 

0.914E-02 

0.424E+00 

0.215E-02 

0.304E-02 

0.728E+00 

0.595E-03 

-0.137E-03 

0.103E+01 

0.557E-04 

-0.265E-04 

0.134E+01 

0.448E-04 

- 

Table  K.4  Reynolds  terms  from  multiple  overheat  data 


y 

Run  009 

Run  010 

0.840E-01 

0.404E-02 

-0.704E-01 

0.301E+00 

0.441E-02 

0.468E-03 

0.517E+00 

0.274E-02 

-0.333E-03 

0.734E+00 

0.187E-02 

-0.438E-02 

0.950E+00 

0.864E-03 

-0.198E-02 
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